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Abstract. Using ab-initio photoionization models of giant H ii regions, we test methods for abundance determinations based on 
a direct measurement of the electron temperature, now possible even for moderate to high-metallicity objects. We find that, for 
metallicities larger than solar, the computed abundances deviate systematically from the real ones, generally by larger amounts 
for more metal-rich H ii regions. We discuss the reasons for this, and present diagrams allowing the reader to better understand 
the various factors coming into play. We briefly discuss less classical methods to derive abundances in metal-rich H ii regions. 
In particular, we comment on the interest of the oxygen and carbon recombination lines. We also show that, contrary to the case 
of giant H ii regions, the physical conditions in bright extragalactic planetary nebulae are such that their chemical composition 
can be accurately derived even at high metallicities. Thus, extragalactic planetary nebulae are promising potential probes of the 
metallicity of the interstellar medium in the internal parts of spiral galaxies as well as in metal-rich elliptical galaxies. 
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1. Introduction 

Giant H ii regions were the first indicators of the presence of 
abundance gradients across the face of galaxies (Searle 1971, 
Shields 1974). Many other abundance indicators have been 
used since then, such as supernova remnants (Dopita et al. 
1980, Smith et al 1993), B-supergiants (Monteverde et al. 1997, 
2000, Urbaneja et al. 2003), and stellar clusters (Ma et al. 
2004, Tiede et al. 2004). A future promising indicator is pro- 
vided by AGB stars (Cioni & Habing 2003). However, H ii re- 
gions remain the most popular tool, both because they are rel- 
atively easy to observe and because the interpretation of their 
spectra in terms of abundances is a priori straightforward. A 
large number of studies based on giant H ii regions has es- 
tablished the existence of abundance gradients in spiral galax- 
ies (see e.g. Zaritsky, Kennicutt & Huchra 1994 for a review). 
The chemical abundances in low metallicity H ii regions are 
generally considered reliable, provided that the spectroscopic 
data are of good quality, allowing an accurate measurement of 
the electron temperature Te from forbidden line ratios such as 
[O m] /14363/5007. In cases where a direct determination of 
the electron temperature is not possible (because the H ii re- 
gions are faint or metal-rich or because the spectra are not very 
deep), the abundances are estimated from statistical methods 
based on strong lines only. These statistical methods need to 
be calibrated (see e.g. Pagel et al. 1979 for an introduction 
of such methods or Stasiriska 2004 for a recent discussion). 
Many calibrations have been proposed. Discrepancies between 
abundances derived using the same data but different calibra- 
tions can be of a factor around 3 (see e.g. Pindao et al 2002). 
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Calibration is a difficult task which requires a sufficient num- 
ber of representative H ii regions and either detailed modelling 
(taking into account geometrical effects and accounting for all 
the relevant data) or a thorough comparison with independently 
derived, accurate abundances, such as abundances from inter- 
stellar absorption lines (Pilyugin 2003). So far, the last word 
has not been said on these matters. 

With the advent of very large telescopes, it is now possible 
to detect and measure the intensities of very weak lines, includ- 
ing [O III] .14363, [S III] ^6312, or [N n] A5155 which are used 
for electron temperature measurements. It is thus hoped that 
it will be possible to derive accurate abundances in metal-rich 
H II regions using classical Tg-based methods. Indeed, a few 
such measurements have already been published, with the con- 
clusion that the metallicity is not as high as was thought before 
(Castellanos et al. 2002, Kennicutt et al. 2003, Garnett et al. 
2004, Bresolin et al. 2004). 

That the metallicity gradients in the inner parts of disks of 
galaxies might not be as steep as previously thought is also sug- 
gested by other indicators such as infra-red lines in H n regions 
(Willner & Nelson-Patel 2002), B stars (Smartt et al., 2001, 
Munn et al. 2004) and even planetary nebulae (Gorny et al. 
2004). However, the evidence is still based on small data sam- 
ples, or on methods that are not free from criticism. It is thus 
important to assess the reliability of the methods used in the 
case of giant H ii regions and to propose alternatives. It has al- 
ready been shown by Stasiriska (1978a, 2002) that abundances 
derived by Tg-based methods are likely to be biased towards 
low metallicities when applied to high metallicity H n regions. 

This paper develops on this aspect. For this, we compute 
series of photoionization models representative of giant H n 
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regions and analyse them with classical Te-based methods for 
abundance determinations. We explain the reasons for the bias 
that occurs for metal-rich H n regions when these methods are 
apphed. We briefly discuss other direct methods of abundance 
analysis. We next consider extragalactic planetary nebulae as 
abundance indicators and show that in this case Te-based abun- 
dance determinations are devoid of bias. We finally come back 
to the case of giant H n regions, and show that their optical 
spectra are extremely sensitive to local perturbations. The main 
points of this study are summarized in Sect. 6. 

2. Tg-based abundance determinations of model 
giant H ii regions 

2. 1. Photoionization models of giant H n regions 

In order to test the reliability of Tg-based abundance determi- 
nations, we construct series of photoionisation models in which 
we vary the metallicity. The model nebulae are spherical with 
a central ionizing source. Since our paper is purely method- 
ological it is sufiicient to consider ionization by blackbodies. 
We consider luminosities that lead to a total number of ioniz- 
ing photons !2(H°) between 4 x and 4 x 10^^ s"\ which 
roughly corresponds to exciting star clusters in giant H n re- 
gions. We consider three geometries: i) a homogeneous sphere 
of hydrogen density n(H) = 50 cm"^; ii) a geometrically thin 
bubble of density n(H) =100 cm~^; iii) a core-halo geometry 
with an inner region at a constant density n(H) - 1000 cm"^ 
surrounded by a halo where the density decreases as the square 
of the distance to the ionizing source. 

The elemental abundances are varied in lockstep with those 
in the Sun as compiled by Lodders (2003), except for nitrogen 
(we assume that N/O varies proportionaUy to O/H and has the 
solar value when O/H is solar) and helium (assumed constant 
and equal to 0.085). Each series consists of 10 models with 
oxygen abundances going from log O/H -n 12 = 8.3 to 9.2 by 
steps of 0.1 dex. In some series, the abundances of Mg, Si and 
Fe are multiplied by a factor of 0.01 to simulate the efl^ects of 
depletion onto dust grains. 

The photoionization models are constructed with the code 
PHOTO as described in Stasihska & Leitherer (1996) with 
atomic data updated as indicated in the Appendix of the present 
paper. All the models are ionization bounded (except the core- 
halo ones) 

Table 1 Usts the series of models considered in this work. 

2.2. Tiie biases 

We applied to our models the analysis techniques that are used 
for real objects. We derived the electron densities from the [S ii] 
/16731/6717 line ratio and the temperatures Tr[0 iii], rr[N ii] 
and T,[0 n], respectively, from the [O ni] /i4363/5007, [N nj 
A5755/6584 and [O ii] i3726,3729/7320,7330 ratios. Then we 
derived the oxygen abundance using standard Tg-based tech- 
niques and assuming that O/H = (0+ + 0^+)/H"^. Note that, as 
stressed by e.g. Stasihska (1978a), these temperature estimates 
are not identical to the mean ionic temperatures for these ions 
(defined as r(0++) = / re«(0++)«edV/ /«(0++)nedyforO-'-', 



and similarly for other ions). The differences may be signifi- 
cant in the case of large temperature gradients or fluctuations. 
In that case, the temperature characteristic of an emission line, 
e.g. ri[0 m] for the [O iii] /15007 line (whose definition is given 
in Stasijiska 1978a) is difl'erent from both the line ratio temper- 
ature Tr[0 m] and the mean ionic temperature 7'(0'^"^). 

The atomic data used in the abundance derivations are ex- 
actly the same as those used in the model computation. This 
allows us to investigate the biases due to the method, indepen- 
dent of uncertainties in the atomic data (which are expected to 
be small, of the order of 10% or less for the ions we consider). 

Fig. 1 shows the derived values of log O/H +12 as a func- 
tion of the input values, for all the models of giant H n regions 
considered. In panel a, the abundances were obtained assum- 
ing that ri[0 ni] = TAO iii] and ri[0 ii] = TA'N ii], i.e. what an 
observer would do if he could measure the [O iii] /14363 and 
[N n] /t5755 line intensities. In panel b, it was assumed that 
Ti[N II] = TAN n] and that ri[0 iii] = (ri[N ii] - 3000) / 0.7, a 
formula commonly used (e.g. Campbell et al. 1986, Kennicutt 
et al. 2003, Gamett et al. 2004a, Bresolin et al. 2004), based 
on the grid of photoionization models by Stasihska (1982). It 
is seen that, as long as the metallicity is low, the derived O/H 
value is very close to the input one (the sUght deviations seen 
in Fig. 2a at low metallicity are due to the fact that the adopted 
empirical relation between ri[0 iii] and ^^[N u] is fulfilled by 
the models only in first approximation and with some scatter). 
Important deviations appear around log O/H -i- 12 = 8.6, and 
may become huge as the metallicity increases. In the case of 
Fig. la, all the derived values are smaller than the input ones, 
sometimes by enormous factors. This means that, if such metal- 
rich H II regions exist, and if observations allow one to measure 
the corresponding line ratios, the method applied in Fig. la will 
always leads to sub-solar derived oxygen abundances! On the 
other hand, the method applied in Fig. lb will either underes- 
timate or overestimate the oxygen abundance. The bias may 
exceed 0.2 dex when the true abundance is larger than the so- 
lar one (log O/H + 12 = 8.7, Allende Prieto et al. 2001). In 
the most extreme cases, however, the transauroral lines are so 
weak that they might not be observable even with very high 
signal-to-noise spectroscopy. In Fig. la we represent by smaller 
symbols the models for which [O m] A4363flifi is smaller than 
10 while in Fig. lb smaller symbols correspond to models 
for which [N n] A5155/H/3 is smaller than 10 This is the flux 
limit attained by Esteban et al. 2002 with 3 hours of integration 
on a 4 m telscope for the brightest H n region in M33. We thus 
see that the best observations available today may lead to con- 
siderably biased abundance determinations. In extreme cases, 
one may use the measured T^O iii] value to judge the like- 
liness of the inferred O/H value. For example, if we find log 
O/H + 12 = 8.3 with the method used in Fig. la while ^^[O m] 
is found to be equal to 5000 - 7000 K, the true abundance is 
likely larger than solar, as can be inferred from Fig. 2a, to be 
described later. However, less extreme cases are more diflicult 
to judge, and would require proper photoionization modelUng 
(but see Sect. 4.4). 
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Table 1. Definition of the model sequences 



name 


^eff 


geometry 


n(H) 


e(H") 


f/in 






depletion 


colour 


symbol 




K 




cm"'' 


photons s"' 














45S1 


45000 


filled sphere 


50 




3.510' 


3.510- 


-3 


no 


red 


squares 


45S1D 


45000 


filled sphere 


50 


410« 


3.510' 


3.510- 


-3 


0.01 


red 


curved squares 


45S2 


45000 


filled sphere 


50 


410'' 


1.610^ 


1.610- 


-2 


no 


red 


circles 


45B 


45000 


bubble 


100 


410"^' 


2.310-3 


2.210- 


-3 


no 


red 


triangles 


45BD 


45000 


bubble 


100 


410"" 


2.310-3 


2.210- 


-3 


0.01 


red 


diamonds 


45CH 


45000 


core-halo 


1000-50 


410^' 


4.610^ 


2.010- 


-1 


no 


red 


stars 


50S1 


50000 


filled sphere 


50 


41049 


4.210' 


4.210- 


-3 


no 


black 


squares 


SOS ID 


50000 


filled sphere 


50 


410« 


4.210' 


4.210- 


-3 


0.01 


black 


curved squares 


50BD 


50000 


bubble 


100 


410^^' 


2.710-3 


2.710- 


-3 


0.01 


black 


diamonds 


lOOPN 


100000 


planetary nebula 


1000 


110*' 


1.610-- 


1.810- 


-3 


no 


blue 


filled circles 


150PN 


150000 


planetary nebula 


10000 


510*^ 


5.5 10-- 


4.810- 


-3 


no 


blue 


filled squares 
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Fig. 1. The oxygen abundance derived by Te-based methods for sequences of model H ii regions as a function of the oxygen abun- 
dance in the models. Each sequence is represented by a different symbol, indicated in Table 1 . In panel a, O/H is computed using 
7'i[0 III] and 7'i[N ii]. Smaller symbols correspond to models with [O iii] /14363/H;6 < IQ-"*. In panel b, O/H is computed using 
Ti-[N n] only. Smaller symbols correspond to models with [N ii] A5155/H/3 < IQ-"*. This figure shows that, at high metallicities, 
the derived O/H values deviate strongly from the true ones. 



3. The cause of the abundance biases at high 
metallicities 

3.1. Temperature gradients inside the nebulae 

It has long been known (Stasinska 1978b, 1980, Garnett 1992) 
that at high metallicities one expects large temperature gradi- 
ents inside H ii regions, due to the fact that the inner zone con- 
tains more efficient coolants (principally O^^ through its infra- 
red lines) than the outer zones. In the following, we elaborate 
on the causes and consequences of these temperature gradients. 

Fig. 2a shows ^[ [O iii] as a function of log 0/H H-12 while 
Fig. 2b shows T, [N ii] as a function of log 0/H -1-12. We see that 
both temperatures decrease with increasing abundances. This is 
due to the enhanced cooling by heavy elements. At high metal- 
licities, however, neither rr[0 m] nor ^[ [N ii] are close to the 



corresponding mean ionic temperatures, as seen in Figs. 3a and 
3b which show T,[0 iii]/r(0++) and r,-[N ii]/r(N+), respec- 
tively, as a function of log O/H -1-12. The values of ^^[N ii] and 
rr[0 m] strongly overestimate the values of T{N^) and T{0^^) 
at metallicities above log O/H H-12 =8.8. The models show that 
they also significantly overestimate the values of 7'i[N n] and 
ri[Ora]. 

It is useful to visualize why temperature gradients are pre- 
dicted to be so large in high metallicity H ii regions. Let us take 
as an example model 45S2 at solar metallicity, and consider one 
point in the high excitation region defined by a local O^^/O^ 
ratio of 0.9, and one point in the low excitation region, defined 
by a local O^^/O^ ratio of 0. 1 . We plot in Fig. 4a the Tg depen- 
dence of log r and of log A in the high excitation region (F and 
A are, respectively, the total energy gains and losses by the gas. 
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Fig. 2. Panel a: T^[0 iii] as a function of O/H for our model H n regions. Panel b: T^[N ii] as a function of O/H. Same symbols as 
in Fig. 1. 
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Fig. 3. Panel a: rr[0 m]/r(0^^) as a function of O/H for our model H ii regions. Panel b: Tr[N ii]/r(N^) as a function of O/H. 
Same symbols as in Fig. 1. This figure illustrates that, at high metallicity, r, [0 iii] and Tr[N ii] are not representative of the O^^ 
and zones, respectively. 



per unit time, per electron and per hydrogen particle). Both F 
and A are in ergs cm-* s"' . Fig. 4b shows the same but for the 
low excitation zone, defined by a local O^^/O^ ratio of 0.1. 
In Fig. 4a, we also plot the emission in the O^^ lines per unit 
time, per electron and per hydrogen atom. The various lines are 
represented in different colours (visible in the on-line version), 
with the following coding: green: O iii] /11661, blue: [O iii] 
.14363, cyan: [O iii] /14949, 5007, magenta: [Oiii] ^52jum, yel- 
low: [Om] /I88yum. In this figure, all the curves include the 
change in the ionization balance due to the variation in electron 



temperature. We see that the oxygen lines are a major contrib- 
utor to the cooling in the high excitation region. The electron 
temperature at each point in a nebula is defined by the con- 
dition F = A. It can easily be shown (see Osterbrock 1989 or 
Stasinska 2002) that F is roughly proportional to the effective 
temperature of the exciting stars and to the recombination coef- 
ficient of hydrogen, which varies approximately like ' . The 
temperature dependence of A results from the sum of the con- 
tribution of the diff'erent cooling lines. At high temperatures, 
the losses are dominated by [O iii] /15007 and increase with 
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Te [lO^K] Te [lO^K] 

Fig. 4. Energy losses (thick increasing curve) and energy gains (thick decreasing curve) as a function of T^, in the 45S2 H ii 
region model at solar metallicity. Panel a: high ionization zone (O^^/O^ - 0.9). The thin curves represent the emissivities O m] 
/11661 (green), [O iii] ^4363 (blue), [O m] /14949, 5007 (cyan), [Om] A52^im (magenta), [Om] AUjim (yellow). Panel b: low 
ionization zone (O^^/O^ - 0.1). The thin curves represent the emissivities of [O ii] /13729 (green), [O ii] /13726 (blue), [N ii] 
A\2l jjm (cyan), [N ii] A205 fim (magenta), [Si n] A35 ;um (yellow). 



Te, while at low temperatures they are dominated by the in- 
frared oxygen lines, and are almost independent of T^. In the 
low excitation region, oxygen is less dominant in the loss pro- 
cesses. The thin curves in Fig. 4b represent the emission in the 
most important cooling lines, with the following colour coding: 
green: [O ii] ^3729, blue: [O ii] ^3726, cyan: [N ii] /1121 jjm, 
magenta: [N ii] A2Q5 fim, yellow: [Si ii] A35 fim. Note that the 
temperature in the low ionization zone, 8700 K, is higher than 
in the high ionization zone, 5500 K, both because cooling is 
less efficient and because heating is more efficient due to the 
hardening of the radiation field in the outskirts of the nebula. 
At low temperatures, [Si n] A35 fim heavily contributes to the 
losses in the low ionization zone, which implies that the tem- 
perature of the gas is strongly dependent on the degree of de- 
pletion of Si onto grains (as already pointed out by Henry 1993 
and Shields & Kennicutt 1995). For a higher metallicity than 
represented in Fig. 4, the entire loss curves will be shifted up- 
wards by a quantity roughly equal to the log of the metallicity 
increase, and the electron temperature, defined by the condi- 
tion F = A, will decrease. For an increase in metallicity by a 
factor of two from the solar value, the temperature will drop 
from Te = 5500 K to - 1500 K in the high ionization zone, 
and from - 8700 K to - 5000 K in the low ionization 
zone, resulting in a larger gradient than in the case of solar 
metallicity. Note that, at the equilibrium temperature for twice 
solar metallicity, the optical and ultraviolet lines from O^^ ions 



are practically not emitted in the high ionization region. Most 
of the emission in these lines actually comes from a layer con- 
taining mostly O^, so that the temperature is sufficiently large 
to allow collisional excitation of these lines. Therefore ^^[O iii] 
and ri[0 III] are not representative of the O^^ zone. 



Fig. 4 also allows one to understand why, in the case of 
large temperature gradients, Tj{0 iii] is larger than T\[0 iii], the 
temperature characteristic of the [O iii] A5QQ1 line: the [O m] 
/14363/5007 ratio is largely weighted by the zones of highest 
temperatures. 



The expected low temperatures and large temperature 
gradients at high metallicities have dramatic consequences 
for abundance determinations from optical lines. The de- 
rived abundances can be under- or overestimated according to 
whether the temperatures adopted to derive the line emissivi- 
ties actually over- or underestimate the true line temperatures. 
The second case may occur if the temperature in the low ex- 
citation zone is estimated not directly from observations but 
from inappropriate model grids, as is the case in Fig. lb. When 
both Ti [O III] and T'i [N ii] can be measured, the -method will 
strongly underestimate the abundances (as seen in Fig. la). 
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^ 5 10 15 20 
Te [lO^K] 

Fig. 5. The Te dependence of the [N ii] /15755/6584 ratio (thin 
curve) and of the [O ii] /l3726,3729/7320,7330 ratio (thick 
curve). The monotonic curves ignore the contribution of re- 
combination. The non-monotonic curves have been computed 
taking into account recombination from N"^"^ and O"^"^ and as- 
suming that these ions are present in the same proportion as N"^ 
and 0+. 

3.2. Contribution of recombination to collisionally 
excited lines. 

Another factor that affects abundance determinations is the 
contribution of recombination to the emissivity of forbidden 
hnes. This has been shown to be important in planetary nebulae 
with cool, hydrogen-poor inclusions (Liu et al. 2000, Tsamis et 
al. 2003a). The importance of this effect in metal-rich H n re- 
gions has not been fully recognized, although it has been con- 
sidered by Kennicutt et al. 2003. These authors have corrected 
the [On] /17320,7330 line emissivity for recombination using 
the formula given by Liu et al. (2000): [O n] ^7320,7330/^8 = 
9.36 (T^/loy-'^ 0++/H+ (which nb. is vahd only for T^/IO"^ 
between 0.5 and 1). For high metallicity H n regions this 
formula may lead to considerable errors. First, if O^'^/H^ is 
obtained from colUsionally excited hnes, it will be overesti- 
mated. Second, the electron temperature characteristic of the 
Hy6 emission is different from the one characteristic of the 
[O nj i7320,7330 emission. 

In Fig. 5 we illustrate the effect of recombinations 
on the temperatures derived from [N n] /15755/6584 and 
[O II] ^3726,3729/7320,7330. The thin curves represent the 
[N II] /15755/6584 ratio and the thick ones the [O n] 
/l3726,3729/7320,7330 ratio. The non-monotonic curves have 
been computed taking into account recombination from N"^"^ 
and O"^"^ and assuming that these ions are present in the same 
proportion as and O'*". The monotonic curves have been 
computed ignoring recombination (these are the ones gener- 



ally used for temperature determinations). We see that both ra- 
tios can be strongly affected by recombination. Below a certain 
temperature, these ratios are even dominated by recombination. 
The minimum of the [O n] /l3726,3729/7320,7330 ratio occurs 
at an electron temperature of around 6000 K, which can be ob- 
tained in the zones for metallicities of the order of log O/H 
-H 12 = 8.75 - 8.9. For [N ii] i5755/6584, the minimum occurs 
at 5000 K, corresponding to O/H + 12 of 8.9 - 9. 

Note that forbidden line from and S^^ are also expected 
to be affected by recombination. However, atomic data that 
would make it possible to compute the contribution of recom- 
bination are not yet available for these ions. 

4. How to eliminate the biases in abundance 
determinations 

4.1. Derive abundances from far infrared lines 

Infrared Une emissivities of metals are almost independent of 
temperature, and thus in principle more suitable than optical 
hnes to derive abundances in the case of metal-rich H n regions. 
Their interpretation poses other problems, though (Simpson et 
al. 2004). So far, studies of extragalactic H n regions using far 
infrared lines are scarce. Willner & Nelson-Patel (2002) have 
used the measurements of [Ne n] and [Ne m] infrared lines in 
25 H II regions in M33 to derive the neon abundance gradient. 
They found that this gradient is much shallower than the oxy- 
gen abundance gradient derived from optical data. However, 
their Ne/H ratios were obtained assuming a temperature of 
10"* K. This assumption creates a bias, since high metallicity 
H n regions will have smaller than low metalhcity ones. 
In the case when the average electron temperature is 5000 K, 
which corresponds to our H ii region models having log O/H 
-1-12 between 8.7 and 8.9, the procedure adopted by Willner & 
Nelson-Patel underestimates the neon abundance by a factor 2. 
If the average temperature of the H ii gas could be measured 
accurately then empirical methods using infrared hnes would 
provide accurate neon abundances in H n regions. 

4.2. Measure the average electron temperature of the 
entire nebula 

The Balmer jump is sensitive to electron temperature and is 
very large at low electron temperatures. Therefore, in principle, 
it should be a good temperature indicator at high metallicities. 
Fig. 6 shows the value of the Balmer jump, defined as the ratio 
of nebular fluxes at 3642 and 3648 A, for our various sequences 
of models. The measurement of the Balmer jump in suppos- 
edly metal-rich extragalactic H n regions should be an impor- 
tant element in deriving accurate abundances. Unfortunately, 
this measurement is difficult, and the contribution of the stellar 
populations must be adequately subtracted. 

4.3. Derive abundances from pure recombination lines 

At metallicities above solar, our models of giant H n regions 
predict for the recombination hnes C n /14267 and O n A465 1 
intensities of up to one hundredth of Hj8. This is much stronger 
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Fig. 7. The intensities of some weak lines as a function of O/H for our model H ii regions. Panel a: O ii /I4651/Hj6. Panel b: C ii 
i4267/H^. Panel c: [O iii] /14363/H/3. Panel d: [N ii] ^5755/H/3. 



than the predicted intensities for the [O m] /14363 or [N ii] 
/15755 lines, as seen in Fig. 7. The intensity ratios of recom- 
bination lines have the virtue of being only little dependent on 
the electron temperature, therefore abundances obtained from 
recombination lines are considered accurate. Some bias may 
exist, however, in the case of nebulae with important tempera- 
ture stratification, since the C^^ and O^^ ions, which generate 
these two lines, are present only in the high ionization zones 
while H/3 is emitted everywhere. 

It must be noted, however, that the interpretation of re- 
combination lines of heavy elements involves other problems. 
In planetary nebulae, abundances derived from recombination 
lines are systematically larger than abundances derived from 
collisionally excited lines by factors from 1.5 to 10 (see e.g. 
Liu 2002 for a review). Recombination lines have also been de- 



tected in a few H ii regions (e.g. Peimbert et al. 1993, Esteban 
et al. 1998, 1999, 2002, Tsamis et al. 2003b, Gai-cia-Rojas 
et al. 2004). Here again they lead to higher abundances than 
collisionally excited lines, although not to the same extreme 
as in planetary nebulae. Temperature stratification cannot be 
the explanation, since this discrepancy has also been found in 
Magellanic Clouds H ii regions that are not expected to harbour 
strong temperature gradients because of their subsolar metal- 
licities. Tsamis et al. (2003b) make the point that this discrep- 
ancy cannot be attributed to temperature fluctuations such as 
described by Peimbert (1967). Presently, the most favoured ex- 
planation is chemical composition inhomogeneities (e.g. Liu et 
al. 2000, Tsamis et al. 2003b) involving extremely metal-rich 
clumps, whose nature is not necessarily the same in planetary 
nebulae and in H ii regions. 
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Fig. 6. The Balmer jump as a function of O/H for our model 
H n regions. 



Therefore recombination Hnes, while offering important 
additional clues to the heavy element abundances in H ii re- 
gions, are not easy to interpret. 

4.4. Tailored photoionization modelling 

In principle, photoionization models take into account all the 
processes that affect the thermal balance at each point in the 
nebula. They also take into account recombination for the emis- 
sion of C, N, O and Ne, at least PHOTO and CLOUDY (Ferland 
1996) do. Therefore, if one fits the observed emission line ra- 
tios with a tailored model, the concerns expressed above are 
not relevant. To simplify the problem, let us assume that only 
oxygen plays an important role in the thermal balance. The 
main parameters defining a model are the mean effective tem- 
perature of the stars, the characteristic ionization parameter, 
and the oxygen abundance, i.e. 3 parameters. The main ob- 
servational constraints ai'e [O iii] /i5007/Hj6, [O iii] /14363/Hj6, 
[O II] ^3727/Hj8 and [O n] A1325/Bfi, i.e. in principle more 
than enough in this simple representation of H ii regions (since 
the measurement of [O iii] /14363 alleviates the double-value 
problem encountered by strong line methods). 

Castellanos et al. (2002) have produced simple tailored 
models for a dozen supposedly metal-rich giant H ii regions 
assuming a simple bubble geometry, and a reasonable solution 
could be found in many cases (although the optimum number 
of constraints was not always available), with the result that the 
metallicity appears to be below solar in the majority. Actually, 
even without detailed photoionization modelling, the fact that 
most of the objects studied by Castellanos et al. (2002) have 
Tr[0 III] and/or T,[N ii] larger than 9000 K imphes that their 
metallicities are not above solar. This can be seen directly on 
our Figs. 2a and 2b (the use of a more reahstic description of 



the stellar radiation field in these figures would not affect this 
claim severely ). Of course, an accurate abundance derivation 
requires a correct model-fitting. 

However, the fact that models are not always able to re- 
produce the optical spectra even in low metallicity giant H n 
regions (e.g. Luridiana et al. 1999, Stasiriska & Schaerer 1999) 
should warn us against being too confident in models for high 
metallicity H ii regions. For the moment, the thermal balance 
of H II regions is not fully understood, and this lack of under- 
standing might have dramatic consequences for the determina- 
tion of abundances in metal-rich H ii regions, where the optical 
forbidden line emissivities depend so strongly on electron tem- 
perature. 

Other factors, such as the density distribution of the nebular 
gas and especially the degree of depletion of metals onto grains 
and the dust content (Henry 1993, Shields & Kennicutt 1995), 
also affect the spectrum of H ii regions and call for a larger 
number of observational constraints (e.g. mid- and far-infrared 
data). 

If indeed large temperature gradients in metal rich H n re- 
gions exist, the spatial distribution of emission line ratios would 
be a useful constraint, as shown by Stasihska (1980). Because 
of the extremely low temperatures in the high ionization zones, 
collisionally excited lines are expected to be strong only in low 
ionization zones (see Fig. 4). Therefore, contrary to what is 
observed at subsolar metallicities, the [O m] /15007 and He i 
/15876 emitting zones should not be coextensive. Note that this 
property does not depend on the details of the nebular geome- 
try. 

Another important constraint on the models comes from in- 
frared lines. Thanks to ISO and the Spitzer Space Telescope, 
these are becoming available for extragalactic H ii. So far, the 
only tailored modelling of an extragalactic H ii region that uses 
optical spectra and images as well as infrared line intensities 
are the ones developed by Garnett et al. (2004b) for a suppos- 
edly metal-rich giant H ii region in M 51 and by Jamet et al. 
(2004) for a metal-poor giant H ii region in M33. In the first 
case, the optical data alone did not provide any direct infor- 
mation on the electron temperature. Constructing photoioniza- 
tion models, the authors estimated the oxygen abundance that 
was able to reproduce the observed [O iii] /l5007/[Ora] /IBS jum 
ratio. However, the observational constraints for this H ii re- 
gion were not many and the gas density distribution was not 
explored. On the other hand, Jamet et al. (2004) studied an H n 
region with measured [O iii] /l4363/[0 iii] /15007 ratio, and con- 
structed photoionization models constrained by the observed 
population of the exciting stars and by the observed nebular 
geometry. They were not able to reproduce concomitantly the 
optical and infrared line fluxes, implying that the uncertainty 
in abundances is larger than thought, even in the case of low 
metallicities. 

5. The use of planetary nebulae as substitutes for 
H II regions 

As we have seen above, the reason why, in truly metal-rich H n 
regions, abundances derived using empirical Tg-based methods 
are likely to be biased is due to the very low electron temper- 
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ature expected in a large fraction of the volume of the H ii re- 
gion. At such a low temperature, collisional excitation of opti- 
cal lines does not occur, and the integrated emission of [O iii] 
lines is strongly weighted by the hottest zones of the H ii re- 
gion, resulting in the [O m] lines being actually emitted where 
O^^ is present only in small amounts. 

Such a low temperature in the O^^ zone of metal-rich nebu- 
lae would not occur if the nebulae were ionized by much hotter 
stars (increased heating) and if they were denser (reduced col- 
lisional deexcitation of infrared cooling lines). Extragalactic 
planetary nebulae that are bright enough for high signal-to- 
noise spectroscopy present both these advantages. It can be 
shown, by using simple photoionization models for expand- 
ing planetary nebulae surrounding evolving post-AGB stars 
(Stasinska et al. 2004) that the most luminous planetary neb- 
ulae (i.e. the ones for which the largest number of important 
diagnostic lines can be observed) are nebulae that have central 
stars with temperatures above 1 00000 K, and hydrogen densi- 
ties above 1000 cm"-'. 

We have constructed two series of models of planetary neb- 
ulae fulfilling these requirements. They are presented at the 
bottom of Table 1 . The set of elemental abundances is identical 
to the one adopted for the giant H ii regions models. 

Figs. 8a and 8b show the energy gains and losses in the 
same way as Fig. 4a and 4b but for the solar abundance plane- 
tary nebula model in the 150PN series instead of the model H ii 
region. These figures show that the temperatures in the O^^ 
and zones are similar, and that the [O iii] A5QQ1 and [O iii] 
/14363 lines are easily emitted in the O^^ zone, so that Tr[0 iii] 
is really representative of the temperature in the O^^ zone (and 
also in the zone to a good approximation). Comparing with 
Fig. 4, r has increased by about a factor 3 and A has decreased 
at the low temperature end. These changes lead to an equilib- 
rium temperature of 11500K (instead of 5500K for the giant 
H n region model at same metallicity shown in Fig. 4). At 
11500 K, the [O m] ^5007 and [O iii] .14363 ai-e emitted effi- 
ciently. 

Despite the orders of magnitude difference between intrin- 
sic luminosites of planetary nebulae and giant H ii regions, the 
[O III] /14363 line is actually stronger in bright planetary neb- 
ulae than in giant H n regions of the same metallicity when 
metallicities are larger than solar. This is shown in Fig. 9, where 
we plot the [O iii] /14363 fluxes from our model giant H ii re- 
gions (panel a) and planetary nebulae (panel b), supposed lo- 
cated at a distance of 1 Mpc. 

Similarly to what we did with our H ii region models, 
we derived the abundances from our planetary nebula mod- 
els using the classical Te-based method. Here, we assumed 
that Ti[0 III] and 7'i[0 ii] are equal to T,^[0 iii]. For oxygen, 
we employed the ionization correction factor commonly used 
for planetary nebulae, i.e. : O/H = (0+ + 0++)/H+ x [(He+ -i- 
He++)/He+]2/3 (Kingsburgh & Barlow 1994). Fig 10 shows the 
derived values of log O/H +12 as a function of the input val- 
ues, for both our series of models of planetary nebulae. It is 
seen that the derived values of log O/H +12 are very close to 
the input ones in the entire metallicity range. For the model 
with higher Teff and higer ionization parameter (model 150PN 
in Table 1), the derived values are systematically lower than the 
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Fig. 10. The oxygen abundance derived by Tg-based methods 
for our model planetary nebulae as a function of the oxygen 
abundance with which the models have been computed. The 
abundances are computed using 7', [0 iii] and 7',[N ii]. 



input ones by 0. 1 dex, because of the inadequacy of the adopted 
ionization correction factor for this precise model. This error is 
however much smaller than the errors in derived abundances 
for giant H ii regions, and occurs only at very high effective 
temperatures and ionization parameters. 

The reason why abundances derived in planetary nebulae 
are so accurate compared to abundances derived in giant H ii 
regions is the larger electron temperature and the smaller in- 
ternal temperature gradients. Note that, given the large values 
of Tg even for the highest metallicity models, the [O ii] /13727 
line is not strongly affected by recombination, contrary to what 
occurs in metal-rich H ii regions. 

There is some concern that, in planetary nebulae, oxygen 
might not trace the abundance of the interstellar medium ex- 
actly, principally due to nuclear and mixing processes in the 
progenitor star (Dopita et al. 1997). If this were true, plane- 
tary nebulae would not be convenient objects for studying the 
metallicity gradients in galaxies. On the other hand, there is no 
such problem with neon (Marigo et al. 2003). It is therefore in- 
teresting to derive the neon abundance for our model planetary 
nebulae. Fig. 1 1 shows the derived values of log Ne/H + 12 as 
a function of the input ones, assuming the commonly used re- 
lation Ne/O - Ne++/0++. We see that the derived Ne/H values 
are within 0.15 dex of the input ones, and that no systematic 
deviation appears with metallicity. For comparison, we show 
in Fig. 12a and b the derived values of log Ne/H -1-12 for our 
model H ii regions adopting the same laws for the temperature 
as for O/H in Figs la and b, respectively. We see that Ne/H de- 
viations from the input values are smaller than in the case of 
O/H, but that they can still be significant. At the highest metal- 
hcities, there is a steep rise in the derived Ne/H, even if O/H 
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Fig. 8. Energy losses (thick increasing curve) and energy gains (thick decreasing curve) as a function of T^, in the 150PN plan- 
etary nebula model at solar metallicity. Panel a: high ionization zone (O^^/O^ - 0.9). The thin curves represent the emissivities 
O m] A\66\ (green), [O m] /14363 (blue), [O m] /14949, 5007 (cyan), [Om] ASl^im (magenta), [Om] AU^im (yellow). Panel b: 
low ionization zone {O^^jO"^ = 0. 1). The thin curves represent the emissivities of [O ii] /13729 (green), [O ii] /13726 (blue), [N ii] 
A\2l yum (cyan), [N ii] A205 fim (magenta), [Si ii] A35 ;um (yellow). 
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Fig. 9. The total flux in the [O iii] /14363 line at a distance of 1 Mpc. Panel a: for our model H ii regions. Panel b: for our model 
planetary nebulae. 



is underestimated (compare Figs la,b and Figs.l2a,b). This is dominates the production of the [Ne iii] /13869 line at temper- 
due to photoionization of Ne^ to excited levels of Ne^^, which atures below, say 3000 - 4000 K (of course, in such cases, the 
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Fig. 12. The Ne/H ratio derived by Te-based methods for our model H ii regions as a function of input Ne/H values in our models. 
In panel a, the abundances are computed using T^IO m] and ^[ [N n]. Smaller symbols correspond to models with [O m] A4363/H/3 
< 10 In panel b, they are computed using ri[N ii] only. Smaller symbols correspond to models with [N ii] A5755fli/3 < 10"^. 
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[Ne III] /13869 line is quite weak, of the order of 10 ^ - 10 ^ of 



5.1. A planetary nebula in a giant H ii region 

Real giant H ii regions are obviously more complex than the 
idealized models we presented here. First of all, their geome- 
tries, as can be judged from high resolution images (Walborn 
et al. 2002, Maiz-Apellaniz et al. 2004) can be quite intricate. 
Besides, the region covered by the spectral slit is likely to con- 
tain emission zones that are not only affected by the ionizing 
photons from massive stars, but may contain supernova rem- 
nants or planetary nebulae. These are not necessarily detected 
individually, but they may have an effect on the observed giant 
H II region spectrum. 

Let us consider again Fig. 4a, to understand qualitatively 
the effect that an intervening planetary nebula may have on 
the spectrum of an H ii region. In an H ii region of half solar 
metallicity, the energy loss curve will drop by about 0.2 dex be- 
low the one represented in Fig. 4a (which corresponds to solar 
metallicity). In the planetary nebula, where the heating is larger 
by a factor about three Te will rise to about 13000K (assum- 
ing that the ionization stage is unchanged). The [O iii] /15007 
emissivity will thus rise by about a factor of two and the [O iii] 
/i4363 emissivity by about a factor of 8. On the other hand, in 
an H n region of twice solar metallicity, the same local rise in 
the energy gains will lead to a local increase in the electron 
temperature from about 2000 K to about 5000 K. This implies 
an increase in the [O iii] A5QQ7 emissivity by over 4 orders of 
magnitude and in the [O iii] /14363 emissivity by an even larger 
amount! Thus, at high metallicities, optical spectra of giant H ii 
regions are extremely sensitive to local variations of the heat- 
ing rate. This is likely to induce important biases when deriving 
their chemical composition. 
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Fig. 13. The giant H ii region 45S2 alone (red circles), the planetary nebula 150PN alone (filled blue circles), and the sum of 
both (black stars). Panel a: 7'r[0 m] as a function of O/H. Panel b: the derived value of O/H as a function of the O/H value in the 
model. 



We illustrate this by considering the superimposition of a 
giant H n region represented by our model 45S2 and of a plan- 
etary nebula represented by our model 150PN. We report in 
Fig. 13a the values of 7'r[0 iii] as a function of log O/H + 12 
computed for the giant H ii region alone (red circles), for the 
planetary nebula alone (filled blue circles), and for the sum of 
both (black stars). We see that, at low metallicity, the value of 
^[ [O III] is not affected by the presence of the planetary neb- 
ula. When log O/H -1-12 exceeds 8.7, rr[0 m] starts increasing 
with O/H. At log O/H -1-12 = 9.2, it reaches the value corre- 
sponding to the planetary nebula alone. As a consequence, at 
a metallicity larger than solar, the bias on O/H would be even 
more dramatic than in the case of a pure giant HII region. This 
assumes that one is able to measure the [O iii] /14363 line in- 
tensity. As in Fig. la, smaller symbols correspond to the cases 
where [O iii] /15007/H/3is smaller than 10 We conclude that 
this "enhanced bias" threatens abundance determinations based 
on the highest quality H ii region spectra obtainable today. 

6. Discussion 

Present-day observations allow one to measure line intensities 
in extragalactic giant H ii regions down to 10^4 of HyS. This 
makes it possible to apply Tg-based methods to derive abun- 
dances even in metal-rich H n regions, while in the past only 
strong line methods could be used. Strong line methods are 
not entirely satisfactory. First of all, they are statistical: they 
assume that all the H n regions obey the same empirical re- 
lation between the ionizing radiation field, the nebular geom- 
etry and its metallicity, which may be true only to a certain 
extent. Besides, the different calibrations proposed in the past 
yield significantly different results for the derived abundances, 
sometimes by factors of 3 or more. This is why the possibility 



of using direct, Te-based methods - believed more precise and 
more accurate - has raised considerable interest. 

In this paper we have examined the biases involved in T^- 
based abundance determinations at high metallicities. Using 
photoionization models of giant H ii regions, we have tested 
these methods by applying to the models the same methods 
as used for real objects. We have found that, for log O/H-i-12 
larger than 8.7 (i.e. larger than the solar value), the computed 
abundances may deviate from the real ones by factors exceed- 
ing 0.2 dex. If both 7'r[0 m] and Ti [N ii] values are used for the 
abundance determinations, the computed O/H values are below 
the true ones, because of strong temperature stratification oc- 
curring in our models. If not all the diagnostics are used but 
guesses are made for some ions, the derived abundances may 
be either underestimated or overestimated, again by significant 
amounts. We have discussed the cause for this, and presented 
diagrams allowing the reader to visualize the various factors 
coming into play. The bottom line of the reason for the fail- 
ure of Tg-based methods in the case of metal-rich H ii regions 
is the fact that optical forbidden lines are significantly emitted 
only above a temperature threshold, and that this threshold lies 
among the temperatures expected in metal-rich H ii regions. 

All our analysis is based on ab initio photoionization mod- 
els that have a simple geometrical structure. In real H ii regions, 
the relative importance of the various factors determining the 
electron temperature in the different zones of the nebulae may 
not be exactly the same. One important problem, not consid- 
ered here, is the effect of dust grains on the thermal balance. 

Because of the temperature dependance of optical forbid- 
den lines and the range of expected values of in high metal- 
licity giant H n regions, any local perturbation of the heating 
rate (shocks, conductive heating) is likely to affect Te-based 
abundances more strongly than in the case of metal-poor H n 
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regions. We have shown that even an intervening high metallic- 
ity planetary nebula may affect Te-based abundance diagnostics 
of high metallicity giant H ii regions. 

We have briefly discussed less classical methods to derive 
abundances in metal-rich H ii regions. In particular, we have 
commented on the interest of the oxygen and carbon recombi- 
nation lines which should be more easily measured in metal- 
rich H n regions than the weak forbidden lines used to estimate 
electron temperatures. 

Detailed photoionization model fitting of presumed metal- 
rich H II regions for which a large number of observational 
constraints can be gathered would be important to better test 
our understanding of the thermal balance of such objects and 
the reliablility of empirical abundance determinations in them. 
We note that, for the moment, even low metallicity giant H ii 
regions are not fully understood from this point of view. 

From the astrophyscal point of view, the main impact of 
abundance determinations in metal-rich H ii regions is on the 
determination of abundance gradients in spiral galaxies. While 
we have shown that strong biases are likely to occur, we have 
not discussed the statistical distribution of these biases. As a 
matter of fact, some insight into the true abundances might be 
gained by considering not only the measured abundance gradi- 
ents (or line intensity gradients) but also their dispersion as a 
function of galactocentric distance. Such an approach requires 
large samples of H ii regions with adequate measurements, 
like those secured recently by Cedres & Cepa (2002) for the 
galaxies NGC 5457 and NGC 4395. It also requires an a priori 
knowledge of the estimated abundance dispersion in the inter- 
stellar medium of galaxies at a given galactocentric distance, 
derived from models of chemical evolution of galaxies. 

We have also shown that, contrary to the case of giant H ii 
regions, the physical conditions in bright extragalactic plane- 
tary nebulae are such that their chemical composition can be 
accurately derived even at high metallicities. Thus, extragalac- 
tic planetary nebulae should be investigated as possible probes 
of the metallicity of the interstellar medium in the internal parts 
of spiral galaxies as well as in metal-rich elUptical galaxies. 

Acknowledgements. Part of this work was done while the author was 
participating in the Guillermo Haro workshop 2004 at the Instituto 
Nacional de Astronomia, Optica y Electronica (Puebla, Mexico) 
which provided a very stimulating and friendly atmosphere. The orga- 
nizers of this workshop are gratefully acknowledged. Thanks also to 
Cathy Ramsbottom and Claire Hudson for giving me access to some 
atomic data. The comments of Enrique Perez and of the referee, Fabio 
Bresolin, are gratefully acknowledged. 

References 

AUende Prieto, C, Lambert, D. L., Asplund, M., 2001, ApJ, 556, L63 
Anderson, H., Ballance, C. R, Badnell, N. R., Summers, H. R, 2002, 

J. Rhys. B: At. Mol. Opt. Phys. 35,1613 
Benjamin R. A., Skillman E. D., Smits D. R, 1999, ApJ, 514, 307 
Berrington K. A., Burke R G., 1981, P&SS, 29, 377 
Berrington K. A., Saraph H. E., Tully, J. A., 1998, A&AS, 129, 161 
Berrington,K.A., 1985, JPhB, 18, L395 

Berrington,K.A., Burke,RG., Dufton,RL., Kingston,A.E., 1985, 
ADNDT, 33, 195 

Bresolin,F., Gamett,D.R., Kennicutt,R.C.,Jr., 2004, astro-ph/0407065 . 



Butler, K., Zeippen, C. J., 1994, A&AS 108, 1 
Campbell,A., Terlevich,R., Melnick,!., 1986, MNRAS, 223, 811 
Castellanos,M., Di'az,A.I., Terlevich,E., 2002, MNRAS, 329, 315 
Cedres B., Cepa J., 2002, A&A, 391, 809 
Cioni,M.-R.L., Habing,H.J., 2003, A&A , 402, 133 
Dopita,M.A., Binette,L., Dodorico,S., Benvenuti,R 1984, ApJ, 276, 
653 

Dopita, M. A., Vassiliadis, E., Wood, P. R., Meatheringham, S. J., 

Harrington, J. P., Bohlin, R. C, Ford, H. C, Stecher, T. R, Maran, 

S. R, 1997, ApJ, 474, 188 
Esteban, C, Peimbert, M., Torres-Peimbert, S., Escalante, V., 1998, 

MNRAS, 295, 401 
Esteban, C, Peimbert, M., Torres-Peimbert, S., Rodriguez, M., 2002, 

ApJ, 581, 241 

Esteban,C., Peimbert,M., Torres-Peimbert,S., Garcfa-Rojas,J., 

Rodn'guez.M., 1999, ApJS, 120, 113 
Ferland, G. J., 1996, Hazy, a Brief Introduction to Cloudy, University 

of Kentucky Department of Physics and Astronomy Internal 

Report 

Froese Fischer 2004a, |http://hf8.vuse.vanderbilt.edu] 

Froese Fischer C, Tachiev G., 2004b, ADNDT, 87,^1 

Gomy, S. K., Stasinska, G., Escudero,A. V., Costa, R. D. D. , 2004, 

astro-ph/0409532 
Galavis, M. E., Mendoza, C, Zeippen, C. I., 1998a, A&AS 131, 499 
Galavis, M. E., Mendoza, C, Zeippen, C. I., 1998b, A&AS 133, 245 
Galavis, M., Mendoza, C, Zeippen, C. J., 1995, A&AS 111, 347 
Galavis, M.E., Mendoza, C, Zeippen, C.J., 1997, A&AS 123, 159 
Garcia-Rojas, J., Esteban, C, Peimbert, M., Rodriguez, M., Ruiz, M. 

T, Peimbert, A., 2004, ApJS, 153, 501 
Gamett, D. R. , 1992, AJ, 103, 1330 

Gamett,D.R., Edmunds, M. G., Henry, R. B. C, Paget, B. E. J., 

Skillman, E. D., 2004b, ApJ 128, 2772 
Garnett,D.R., Kennicutt,R.C.,Jr., Bresolin,F., 2004a, ApJ, 607L, 21 
Griffin D. C, Badnell N. R., 2000, JPhB, 33, 4389 
Griffin D. C, Mitnik D. M., Badnell, N. R., 2001, IPhB, 34, 4401 
Griffin,D.C., Badnell,N.R., Pindzola,M.S., 2000, JPhB, 33, 1013 
Henry, R.B.C., 1993, MNRAS, 261, 306 
Hudson, C. E., Bell, K. L., 2004, MNRAS 348, 1275 
Jamet, L., Stasinska, G., Perez, E., et al. 2004 submitted 
Keenan,F.R, Feibelman,W.A., Berrington,K.A., 1992, ApJ, 389, 443 
Kennicutt,R.C.,Jr., Bresolin,F, Garnett,D.R., 2003, Apl, 591, 801 
Kingsburgh, R. L., Barlow, M. J., 1994, MNRAS, 271, 257 
Liu, X.-W, 2002, Rev Mex A&A, 12, 70 

Liu, X.-W., Storey, P. J., Barlow, M. J., Danziger, I. I., Cohen, M., 

Bryce, M., 2000, MNRAS, 312, 585 
Lodders, K. 2003, ApJ, 591, 1220 

Luridiana, V., Peimbert, M., Leitherer, C, 1999, ApJ, 527, 110 
Ma,J., Zhou,X., Chen,J., 2004, A&A , 413, 563 
Mai'z-Apellaniz J., Perez E., Mas-Hesse J. M., 2004, AJ, 128, 1 196 
Marigo, R, Bernard-Salas, J., Pottasch, S. R., Tielens, A. G. G. M., 

Wesselius, R R., 2003, A&A, 409, 619 
McLaughlin B. M., Bell K. L., 2000, JPhB, 33, 597 
McLaughlin B. M., Bell K. L., 2000, JPhB, 33, 597 
Mendoza, C, Zeippen, C. J., 1982a, MNRAS 202, 981 
Mendoza, C, Zeippen, C. J., 1982b, MNRAS 198, 127 
Mendoza, C, Zeippen, C. J., Storey, R J., 1999, A&AS 135, 159 
Mendoza, C, 1983, lAUS, 103, 143 

Merkehs G., Martinson I., Kisielius R., Vilkas M. J., 1999, PhyS, 59, 
122 

Mitnik D. M., Griffin D. C, Badnell, N. R., 2001, JPhB, 34, 4455 
Monteverde,M.I., Herrero,A., Lennon,D.J, 2000, ApJ, 545, 813 
Monteverde,M.I., Herrero,A., Lennon,D.J., Kudritzki,R.-P., 1997, 
ApJ, 474L.107 



14 



Munn,K.E., Dufton,P.L., Smartt,S.J., Hambly,N.C., 2004, A&A , 419, 
713 

Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae and active 

galactic nuclei (University Science Books) 
Pequignot D., Aldrovandi S. M. V., 1976, A&A, 50, 141 
Pagel, B. E. J., Edmunds, M. G., Blackwell, D. A., Chun, M. S., Smith 

G., MNRAS, 189, 95 
Peimbert, M., 1967, ApJ, 150, 825 

Peimbert M., Storey P J., Torres-Peimbert S., 1993, ApJ, 414, 626 
Pelan J., Berrington K. A., 1995, A&AS, 1 10, 209 
Pilyugin, L. S., 2003, A&A, 399, 1003 

Pindao,M., Schaerer,D., GonzalezDelgado,R.M., Stasiiiska,G., 2002, 
A&A, 394, 443 

Pradhan, AH., 1995, in The Analysis of Emission Lines, 1994, Eds.: 

R. Williams and M. Livio, Cambridge University Press, p. 8 
Ramsbottom C. A., Bell K. L., Keenan, E P, 1998, MNRAS, 293, 233 
Ramsbottom C. A., Bell K. L., Stafford R. P, 1996, ADNDT, 63, 57 
Ramsbottom C. A., Berrington K. A., Hibbert A., Bell K. L., 1994, 
PhyS, 50, 246 

Ramsbottom, C. A., Bell, K. L., Keenan, E P,1997, MNRAS 284, 754 

Saraph, E., Storey, P J., 1996, A&AS, 115, 151 

Searle, L., 1971, ApJ, 168, 327 

Shields, G. A. , 1974, ApJ, 193, 335 

Shields, J.C., Kennicutt, R.C., 1995, ApJ, 454, 807 

Sigut T. A. A., Pradhan A. K., 1995, JPhB, 28, 4879 

Simpson, J. P., Rubin, R. H., Colgan S. W. J., Erickson, E. E, Haas, 

M. R., 2004, ApJ, 611, 338 
Smartt,S.J., Venn,K.A., Dufton,PL., Lennon,D.J., Rolleston,W.R.J., 

Keenan,EP, 2001, A&A , 367, 86 
Smith,R.C., Kirshner,RobertP., Blair, WilliamP, Long,KnoxS., 

Winkler,PE 1993, ApJ, 407, 564 
Stafford R. P, Bell K. L., Hibbert A., Wijesundera W. R, 1994, 

MNRAS, 268, 816 
Stasihska, G., 1978a, A&A, 66, 257 
Stasihska, G., 1978b, A&AS, 32, 429 
Stasihska, G., 1980, A&A, 85, 359 
Stasihska, G., 1982, A&AS, 48, 299 

Stasihska, G., 2002, RevMexAA (Serie de Conferencias), 12, 62 
Stasihska, G., 2004, in Cosmochemistry, The Melting Pot of the 

Elements, Eds. C. Esteban, A. Herrero, R. Garcia Lopez, 

Cambridge University Press, p. 1 15 
Stasihska, G., Leitherer K., 1996, ApJS 107, 661 
Stasihska, G., Schaerer, D., 1999, A&A, 351, 72 
Stasihska, G., Vflchez, J. M., Perez E., Gonzalez Delgado, R., 

Corradi, R. M. L., Mampaso, A., Magrini, L. 2004, ESO confer- 
ence proceedings (in press) 
Storey,P J., Hummer,D.G., 1995, MNRAS, 272, 41 
Storey, PL, Zeippen, C.J., 2000, MNRAS, 312, 813 
Tachiev G., Eroese Eischer C, 1999, JPhB, 32, 5805 
Tachiev G. L, Eroese Eischer C, 2002, A&A, 385, 716 
Tayal, S. S., Gupta, G. P, 1999, ApJ 526, 544 

Tayal, S. S., 2000, ApJ, 530, 1091 

Tiede, G. P, Sarajedini, A., Barker, M. K., 2004, "astro-ph/0403525 1 
Tsamis, Y. G., Barlow, M. L, Liu, X.-W., Danziger, 1. L, Storey, P J., 

2003a, MNRAS, 345, 186 
Tsamis, Y. G., Barlow, M. L, Liu, X.-W., Danziger, 1. L, Storey, P J., 

2003b, MNRAS, 338, 687 
Urbaneja,M.A., Herrero,A., Kudritzki,R.P., Bresolin,E., Corral,L.J., 

Puls,J., 2002, A&A , 386.1019 
Walborn, N. R., Mafz-Apellaniz J., Barba R. H., 2002, AJ, 124, 1601 
Willner, S. P, Nelson-Patel, K., 2002, ApJ, 568, 679 
Wilson, N. J., Bell, K. L. , 2002, MNRAS, 337, 1027 
Zaritsky, D., Kennicutt, R.C., Huchra, LP, 1994, ApJ, 420, 87 
Zatsarinny O., Tayal S. S., 2003, ApJS, 148, 575 
Zhang H. L., Graziani M., Pradhan, A. K., 1994, A&A, 283, 319 
Zhang H. L., Sampson D. H., 1993, PhRvA, 47, 208 



15 



Appendix A: Atomic data 

In the version of PHOTO used in this work, the atomic data to 
compute line emissivities have been updated with respect to the 
ones listed in Stasihska & Leitherer (1996). The recombination 
line emissivities for HI are from Storey & Hummer (1995) with 
collisional excitation contribution using data from Anderson et 
al. (2002). The emissivities of Hel lines are from Benjamin et 
al. (1999) and those of Hell lines are from Storey & Hummer 
(1995). For the remaining elements, the references are given in 
Table A. 1 . In this table, the reference Froese Fischer & Tachiev 
(2004a) refers to the internet site http://hf8.vuse.vanderbilt.edu 
After our compilation of atomic data was performed, the paper 
Froese Fischer & Tachiev (2004b) was published, giving tran- 
sition probabilities of all ions of interest. 



Table A.l. Atomic data for line emission 



ion 


collision strengths 


transition probabilities 


C I 


Pequignot & Aldrovandi (1976) 


Galavis et al. ( 1997a), Mendoza & Zeippen (1999) 


CII 


Wilson & Bell (2002) 


Galavis etal. (1998) 


cm 


Berrington (1985, 1992), Burke et al. (1985), 


Tachiev & Froese Fischer (1999) 


CIV 


Griffin et al. (2000) 


Kingston & Hibbert (2001) 


N I 


Berrington & Burke (1981) 


Tachiev & Froese Fischer (2002) 


M TT 


Hudson & Bell (/0U4) 


Galavis et ai. (199 /a), atorey & Zeippen (2U0U) 


NIII 


Stafford etal. (1994) 


Galavis etal. (1998) 


NIV 


Ramsbottom et al. (1994) 


Tachiev & Froese Fischer (1999) 


U I 


Zatsarinny & Tayal (2003) 


Oalavis et al. {lyy /) 


O II 


McLaughlm& Bell (2000) 


Tachiev & Froese Fischer (2002) 


O III 


Aggarwal & Keenan (1999) 


Galavis et al. (1997), Mendoza et al. (1999), Storey & Zeippen (2000) 


yj iV 


Anang ei ai. {lyy'^) 


Lraiavis ei ai. (lyy /d; 


V 


Zhang & Sampson (1993) 


Tachiev & Froese Fischer (1999) 


Nell 


Griffin et al. (2001) 


Froese Fischer & Tachiev (2004a) 


Nelll 


McLaughlin & Bell (2000) 


Galavis et al. (1997a), Storey & Zeippen (2000) 


NelV 


Ramsbottom et al. (1998) 


Merkelis etal. (1999) 


Ne V 


Griffin & Badnell (2000) 


Galavis et al. (1997a) 


Ne VI 


Mitnik etal. (2001) 


Galavis et al. (1998) 


Mgll 


Sigut&Pradhan(1995) 


Froese Fischer & Tachiev (2004a) 


MglV 


Berrington et al. (1998) 




Mg V 


Butler &Zeippen (1994) 


Tachiev & Froese Fischer (2002a) 


Si II 


Pradhan (1995) 


Froese Fischer & Tachiev (2004a) 


SII 


Ramsbottom et al., (1996) 


Froese Fischer & Tachiev (2004a) 


SIII 


Tayal & Gupta (1999) 


Froese Fischer & Tachiev (2004a) 


SIV 


Tayal (2000) 


Froese Fischer & Tachiev (2004a) 


Arll 


Pelan & Berrington (1995) 


Mendoza (1983) 


Arm 


Galavisetal. (1998) 


Mendoza & Zeippen (1982a) 


ArlV 


Ramsbottom et al. (1997) + 2004 private communication 


Mendoza & Zeippen (1982b) 


Ar V 


Galavisetal. (1995) 


Froese Fischer & Tachiev (2004a) 


Ar VI 


Saraph& Storey (1996) 


Froese Fischer & Tachiev (2004a) 



